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Coherent vortical motion has been reported in a wide variety of populations including living 
organisms (bacteria, fishes, human crowds) and synthetic active matter (shaken grains, mixtures 
of biopolymers), yet a unified description of the formation and structure of this pattern remains 
lacking. Here we report the self-organization of motile colloids into a macroscopic steadily rotating 
vortex. Combining physical experiments and numerical simulations, we elucidate this collective be¬ 
havior. We demonstrate that the emergent-vortex structure lives on the verge of a phase separation, 
and single out the very constituents responsible for this state of polar active matter. Building on 
this observation, we establish a continuum theory and lay out a strong foundation for the descrip¬ 
tion of vortical collective motion in a broad class of motile populations constrained by geometrical 
boundaries. 
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Building upon the pioneering work of Vicsek et al. [T], 
physicists, mathematicians and biologists have contem¬ 
plated the self-organization of living-organism groups 
into flocks as an emergent process stemming from simple 
interaction rules at the individual level [IS]. This idea 
has been supported by quantitative trajectory analysis in 
animal groups IMZI, together with a vast number of nu¬ 
merical and theoretical models mil], and more recently 
by the observations of flocking behavior in ensembles of 
non-living motile particles such as shaken grains, active 
colloids, and mixtures of biofilaments and molecular mo¬ 
tors miHiiaiiii. From a physicist’s perspective, these 
various systems are considered as different instances of 
polar active matter, which encompasses any ensemble of 
motile bodies endowed with local velocity-alignment in¬ 
teractions. The current paradigm for flocking physics 
is the following. Active particles are persistent random 
walkers, which when dilute form a homogeneous isotropic 
gas. Upon increasing density, collective motion emerges 
in the form of spatially localized swarms that may cruise 
in a sea of randomly moving particles; further increasing 
density, a homogeneous polar liquid forms and sponta¬ 
neously flows along a well-defined direction miniiisi. 
This picture is the outcome of experiments, simulations 
and theories mostly performed in unbounded or periodic 
domains. 

Beyond this picture, significant attention has been de¬ 
voted over the last five years to confined active mat¬ 
ter j3l 1141426] . Confined active particles have consistently, 
yet not systematically, been reported to self-organize into 
vortex-like structures. However, unlike for our under¬ 
standing of flocking, we are still lacking a unified picture 
to account for the emergence and structure of such vortex 
patterns. This situation is mostly due to the extreme di¬ 
versity in the nature and symmetries of the interactions 
between the active particles that have been hitherto con¬ 
sidered. Do active vortices exist only in finite-size sys¬ 
tems as in the case of bacterial suspensions m, which 
lose this beautiful order and display intermittent turbu¬ 
lent dynamics [27] when unconfined? What are the nec¬ 
essary interactions required to observe and/or engineer 
bona fide stationary swirling states of active matter? 

In this paper, we answer these questions by con¬ 
sidering the impact of geometrical boundaries on the 
collective behavior of motile particles endowed with 
velocity-alignment interactions. Combining quantitative 
experiments on motile colloids, numerical simulations 
and analytical theory, we elucidate the phase behavior of 
polar active matter restrained by geometrical boundaries. 
We use colloidal rollers which, unlike most of the avail¬ 
able biological self-propelled bodies, interact via well 
established dynamical interactions [T]. We first exploit 
this unique model system to show that above a critical 
concentration populations of motile colloids undergo 
a non-equilibrium phase transition from an isotropic 
gaseous state to a novel ordered state where the entire 
population self-organizes into a single heterogeneous 
steadily rotating vortex. This self-organization is not the 


consequence of the finite system size. Rather, this emer¬ 
gent vortex is a genuine state of polar active matter lying 
on the verge of a macroscopic phase separation. This 
novel state is the only ordered phase found when unidi¬ 
rectional directed motion is hindered by convex isotropic 
boundaries. We then demonstrate theoretically that a 
competition between alignment, repulsive interactions 
and confinement is necessary to yield large-scale vortical 
motion in ensembles of motile particles interacting 
via alignment interactions, thereby extending the rel¬ 
evance of our findings to a broad class of active materials. 


RESULTS 

Experiments. The experimental setup is fully de¬ 
scribed in the Methods section and in Figs. I^i and 1^. 
Briefly, we use colloidal rollers powered by the Quincke 
electrorotation mechanism as thoroughly explained in [1] . 
An electric field Eq is applied to insulating colloidal 
beads immersed in a conducting fluid. Above a critical 
field amplitude Eq^ the symmetry of the electric charge 


FIG. 1. Experimental setup, (a) Sketch of the setup. 5- 
micron PMMA colloids roll in a microchannel made of two 
ITO-coated glass slides assembled with double-sided scotch 
tape. An electrokinetic flow confines the rollers at the center 
of the device in a circular chamber of radius Rc. (b) Super¬ 
imposed fluorescence pictures of a dilute ensemble of rollers 
{Eq/Eq = 1.1, 00 = 6 X 10“^). The colloids propel only in¬ 
side a circular disc of radius Rc = 1 nim and follow persistent 
random walks. 
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FIG. 2. Dynamics of an isolated colloidal roller, (a) Local packing fraction 0(r), averaged over the azimuthal angle 0, 
plotted as a function of the radial distance. The dashed line indicates the radius of the circular chamber, (b) Probability 
distribution function of the roller velocities measured from the individual tracking of the trajectories, (c) Autocorrelation of 
the roller velocity (vi(t) • Vi(t + T)) plotted as a function of vqT for packing fractions ranging from 0o = 6 x 10“^ to 0o = 10“^. 
Full line: best exponential ht. (d) Superimposed trajectories of colloidal rollers bouncing off the edge of the conhning circle. 
Time interval: 5.3 ms. {Eq/Eq — 1.1, 0o = 6 x 10“^). Same parameters for the four panels: Rc — 1.4mm, E^jEq — 1.1, 

00 = 6 X 10“^ 


distribution at the bead surface is spontaneously broken. 
As a result, a net electric torque acts on the beads caus¬ 
ing them to rotate at a constant rate around a random 
axis transverse to the electric field [28H3Q] . When the 
colloids sediment, or are electrophoretically driven, onto 
one of the two electrodes, rotation is converted into a 
net rolling motion along a random direction. Here, we 
use PMMA spheres of radius a = 2.4 /rm immersed in a 
hexadecane solution. 

As sketched in Fig. [^, the colloids are handled and 
observed in a microfluidic device made of double-sided 
scotch tape and of two ITO-coated glass slides. The ITO 
layers are used to apply a uniform DC field Eq = Eqz 
in the 2 :-direction, with Eq = 1.6V//im {Eq = I.IEq). 
Importantly, the electric current is nonzero solely in a 
disc-shaped chamber at the center of the main channel. 
As exemplified by the trajectories shown in Fig. and 
in Supplementary Movie i, Quincke rotation is hence 
restrained to this circular region in which the rollers 
are trapped. We henceforth characterize the collective 
dynamics of the roller population for increasing values 
of the colloid packing fraction 0o- 

Individual self-propulsion. For area fractions smaller 
than = 10“^, the ensemble of rollers uniformly ex¬ 
plores the circular confinement as illustrated by the flat 
profile of the local packing fraction averaged along the 
azimuthal direction 0(r) in Fig. [^. The rollers undergo 
uncorrelated persistent random walks as demonstrated in 
Figs. Ht an d^. The probability distribution of the roller 
velocities is isotropic and sharply peaked on the typical 
speed vq = 493±17/im.s“^. In addition, the velocity au¬ 
tocorrelation function decays exponentially at short time 
as expected from a simple model of self-propelled parti¬ 
cles having a constant speed vq and undergoing rotational 
diffusion with a rotational diffusivity D~^ = 0.31 ±0.02 s 
that hardly depends on the area fraction (see Supple¬ 


mentary Diseussion 1). These quantities correspond to a 
persistence length of = vq/D = 160 yum that is about 
a decade smaller than the confinement radius Rc used in 
our experiments: 0.9 mm < Rc < 1.8 mm. 

At long time, due to the collisions on the disc bound¬ 
ary, the velocity autocorrelation function sharply drops 
to 0 as seen in Fig. [^. Unlike swimming cells [26l 1^ . 
self-propelled grains 0 m [22] or autophoretic col¬ 
loids m, dilute ensembles of rollers do not accumulate 
at the boundary. Instead, they bounce off the walls of 
this virtual box as shown in a close-up of a typical roller 
trajectory in Fig. Ht, and in the Supplementary Movie 
1. As a result, the outer region of the circular cham¬ 
ber is depleted, and the local packing fraction vanishes 
as r goes to i?c, Fig. Ht- The repulsion from the edges 
of the circular hole in the microchannel stems from an¬ 
other electrohydrodynamic phenomenon [33]. When an 
electric field is applied, a toroidal flow sketched in Fig. 
Et is osmotically induced by the transport of the electric 
charges at the surface of the insulating adhesive films. 
Consequently, a net inward flow sets in at the vicinity of 
the bottom electrode. As the colloidal rollers are prone to 
reorient in the direction of the local fluid velocity [1] , this 
vortical flow repels the rollers at a distance typically set 
by the channel height H while leaving unchanged the col¬ 
loid trajectories in the center of the disc. This electroki- 
netic flow will be thoroughly characterized elsewhere. 

Collective motion in confinement. As the area 
fraction is increased above 0*, collective motion emerges 
spontaneously at the entire population level. When the 
electric field is applied, large groups of rollers akin to 
the band-shaped swarms reported in [T] form and collide. 
However, unlike what was observed in periodic geome¬ 
tries, the colloidal swarms are merely transient and ulti¬ 
mately self-organize into a single vortex pattern spanning 
the entire confining disc as shown in Fig. Ht and Supple- 
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FIG. 3. Collective-dynamics experiments, (a) Snapshot of a vortex of rollers. The dark dots show the position of one 
half of the ensemble of rollers. The blue vectors represent their instantaneous speed (Re = 1.35 mm, 0o = 5 x 10-^). (b) 
Average polarization plotted versus the average packing fraction for different confinement radii. Open symbols: experiments. 
Full line: best fit from the theory. Filled circles: numerical simulations (b = 3a, Rc = 1mm). (c) Time-averaged polarization 
field {Rc = 1.35 mm, 0o = 5 x 10“^). (d) Time average of the local packing fraction {Rc = 1.35 mm, 0o = 5 x 10“^). (e) 
Time-averaged packing fraction at the center of the disc, normalized by 0* and plotted versus the average packing fraction. 
Error bars: one standard deviation, (f) Fraction of the disc where Ilc^ > 0.5 versus the average packing fraction. Open symbols: 
experiments. Full line: theoretical prediction with no free fitting parameter. Filled circles: numerical simulations {b = 3a, 
Rc = I mm), (g) Radial density profiles plotted as a function of the distance to the disc center r. All the experiments correspond 
to 00 = 0.032 ± 0.002, error bars: la. (h) Open symbols: same data as in (g). The radial density profiles are rescaled by 0* 
and plotted versus the rescaled distance to the center r/Rc. All the profiles are seen to collapse on a single master curve. Filled 
symbols: Numerical simulations. Solid line: theoretical prediction. All the data correspond to EqIEq = 1.1. 


mentary Movie 2. Once formed, the vortex is very robust, 
rotates steadily and retains an axisymmetric shape. In 
order to go beyond this qualitative picture, we measured 
the local colloid velocity field v(r, t) and use it to define 
the polarization field n(r, t) = v/i;o5 which quantifies 
local orientational ordering. The spatial average of 11 
vanishes when a coherent vortex forms, therefore we use 
its projection 11^ = (11 • e^)r,t along the azimuthal di¬ 
rection as a macroscopic order parameter to probe the 
transition from an isotropic gas to a polar-vortex state. 
As illustrated in Fig. |^, n^(0o) displays a sharp bifur¬ 
cation from an isotropic state with FI^ = 0 to a globally 
ordered state with equal probability for left- and right- 
handed vortices above 0o = 0"*". Furthermore, Fig. 
demonstrates that this bifurcation curve does not depend 
on the confinement radius Rc . The vortex pattern is spa¬ 
tially heterogeneous. The order parameter and density 
fields averaged over time are displayed in Figs. Iff and 
p, respectively. At first glance, the system looks phase- 
separated: a dense and ordered polar-liquid ring where all 
the colloids cruise along the azimuthal direction encloses 


a dilute and weakly ordered core at the center of the disc. 
We shall also stress that regardless of the average pack¬ 
ing fraction, the packing fraction in the vortex core is 
measured to be very close to 0"*", the average concentra¬ 
tion below which the population is in a gaseous state, see 
Fig. Iff. This phase-separation picture is consistent with 
the variations of the area occupied by the ordered outer 
ring, Aring, for different confinement radii R^ as shown 
in Fig.l^. We define A^ing as the area of the region where 
the order parameter exceeds 0.5, and none of the results 
reported below depend on this arbitrary choice for the 
definition of the outer-ring region. A^ing also bifurcates 
as 00 exceeds 0*, and increases with Rc. Remarkably, all 
the bifurcation curves collapse on a single master curve 
when Ariiw is rescaled by the overall confinement area 
ttRI^ Fig.j^. In other words, the strongly polarized outer 
ring always occupies the same area fraction irrespective 
of the system size, as would a molecular liquid coexisting 
with a vapor phase at equilibrium. However, if the sys¬ 
tem were genuinely phase-separated, one should be able 
to define an interface between the dense outer ring and 
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the dilute inner core, and this interface should have a con¬ 
stant width. This requirement is not borne out by our 
measurements. The shape of the radial density profiles 
of the rollers in Fig.[^ indeed makes it difficult to unam¬ 
biguously define two homogeneous phases separated by a 
clear interface. Repeating the same experiments in discs 
of increasing radii, we found that the density profiles are 
self-similar. Fig. Si- The width of the region separat¬ 
ing the strongly polarized outer ring from the inner core 
scales with the system size, which is the only character¬ 
istic scale of the vortex patterns. The colloidal vortices 
therefore correspond to a monophasic yet spatially het¬ 
erogeneous liquid state. 

In order to elucidate the physical mechanisms respon¬ 
sible for this intriguing structure, we now introduce a 
theoretical model that we solve both numerically and 
analytically. 

Numerical simulations. The Quincke rollers are elec¬ 
trically powered and move in a viscous fluid, and hence 
interact at a distance both hydro dynamically and electro¬ 
statically. In [1], starting from the Stokes and Maxwell 
equations, we established the equations of motion of a 
dilute ensemble of Quincke rollers within a pairwise ad¬ 
ditive approximation. When isolated, the ith roller lo¬ 
cated at r^ moves at a speed uq along the direction 
Pi = (cos sinopposite to the in-plane component 
of the electrostatic dipole responsible for Quincke rota¬ 
tion [T]. When interacting via contact and electrostatic 
repulsive forces, the roller velocity and orientation are 
related by: 

= t>oPi - y]'Hrep(ri - Tj)- (1) 

Inertia is obviously ignored, and for the sake of simplic¬ 
ity we model all the central forces acting on the colloids 
as an effective hard-disc exclusion Hrep of range b. In 
addition, Oi follows an over damped dynamics in an ef¬ 
fective angular potential capturing both the electrostatic 
and hydrodynamic torques acting on the colloids [1]: 

dtOiit) = de, ^ U{vi - r^; Pi, Pj) + £,i. (2) 

The ^^’s account for rotational diffusion of the rollers. 
They are uncorrelated white noise variables with zero 
mean and variance = 2D6{t — The 

effective potential in Eq. is composed of three terms 
with clear physical interpretations: 

Pi, pj) = A{r) pj ■ pi + B{r) f • p* (3) 
+ C{r)pj ■ (2fr-I) - Pi, 

where r = r/r. The symmetry of these interactions is 
not specific to colloidal rollers and could have been an¬ 
ticipated phenomenologically exploiting both the transla¬ 
tional invariance and the polar symmetry of the surface- 
charge distribution of the colloids m- The first term 


promotes alignment and is such that the effective poten¬ 
tial is minimized when interacting rollers propel along 
the same direction. A{r) is positive, decays exponen¬ 
tially with r/i7, and results both from hydrodynamic and 
electrostatic interactions. The second term gives rise to 
repulsive torques^ and is minimized when the roller ori¬ 
entation points away from its interacting neighbor. B{r) 
also decays exponentially with r/H but solely stems from 
electrostatics. The third term has a less intuitive mean¬ 
ing, and promotes the alignment of along a dipolar 
field oriented along pj . This term is a combination of hy¬ 
drodynamic and electrostatic interactions, and includes 
a long-ranged contribution. 

The functions A{r), B{r), and C{r) are provided in 
the Supplementary Diseussion 2. As it turns out, all the 
physical parameters (roller velocity, field amplitude, fluid 
viscosity, etc.) that are needed to compute their exact 
expressions have been measured, or estimated up to loga¬ 
rithmic corrections, see Supplementary Diseussion 2. We 
are then left with a model having a single free parame¬ 
ter that is the range, 6, of the repulsive forees between 
colloids. We numerically solved this model in circular 
simulation boxes of radius Rc with reflecting boundary 
conditions using an explicit Euler scheme with adaptive 
time-stepping. All the numerical results are discussed 
using the same units as in the experiments to facilitate 
quantitative comparisons. 

The simulations revealed a richer phenomenology than 
the experiments, as captured by the phase diagram in 
Fig. 1^ corresponding to Rc = 0.5 mm. By systemati¬ 
cally varying the range of the repulsive forces and the 
particle concentration, we found that the (0o,^) plane 
is typically divided into three regions. At small pack¬ 
ing fractions, the particles hardly interact and form an 
isotropic gaseous phase. At high fractions, after transient 
dynamics strikingly similar to that observed in the exper¬ 
iments, the rollers self-organize into a macroscopic vortex 
pattern, Eig. and Supplementary Movie 3. However, 
at intermediate densities, we found that collective mo¬ 
tion emerges in the form of a macroscopic swarm cruising 
around the circular box through an ensemble of randomly 
moving particles, Eig. and Supplementary Movie 
These swarms are akin to the band patterns consistently 
reported for polar active particles at the onset of col¬ 
lective motion in periodic domains mm- This seeming 
conflict between our experimental and numerical findings 
is solved by looking at the variations of the swarm length 

with the confinement radius Rc in Eig. [^. We define 

as the correlation length of the density fluctuations 
in the azimuthal direction. The angular extension of 
the swarms ^s/Rc increases linearly with the box radius. 
Therefore, for a given value of the interaction parameters, 
there exists a critical box size above which the popula¬ 
tion undergoes a direct transition from a gaseous to an 
axisymmetric vortex state. Eor b = 3a, which was mea¬ 
sured to be the typical interparticle distance in the polar 
liquid state [1], this critical confinement is Rc = 1mm. 
This value is close to the smallest radius accessible in our 
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FIG. 4. Collective-dynamics simulations, (a) The numerical phase diagram of the conhned population is composed of 
three regions: isotropic gas (low 0o, small b), swarm coexisting with a gaseous phase (intermediate 0o and 5), and vortex state 
(high 00 and b). Rc = 0.5 mm. (b) Snapshot of a vortex state. Numerical simulation for 0o = 0.1, and b = 5a. (c) Snapshot of 
a swarm. Numerical simulation for 0o = 4.5 x 10~^, and b = 2a. (d) Variation of the density correlation length as a function 
of Rc. Above Rc = 1mm, ^ plateaus and a vortex is reached (0o = 3 x 10“^, b = 3a). (e) Four numerical snapshots of rollers 
interacting via: alignment interactions only (A), alignment interactions and repulsive torques (A + B, where the magnitude 
of S is 5 times the experimental value), alignment and excluded volume interactions (A + 5, where the repulsion distance is 
b = 5a), alignment and the C-term in Eq. 3 (A + C). Polarized vortices emerge solely when repulsive couplings exist (A B 
and A + 6). 


experiments where localized swarms were never observed, 
thereby solving the apparent discrepancy with the exper¬ 
imental phenomenology. 

More quantitatively, we systematically compare our 
numerical and experimental measurements in Figs. 
If and for Rc = 1 mm. Even though a number 
of simplifications were needed to establish Eqs. SHI 
and nil], the simulations account very well for the 
sharp bifurcation yielding the vortex patterns as well 
as their self-similar structure. This last point is proven 
quantitatively in Fig. [^, which demonstrates that the 
concentration increases away from the vortex core, 
where 0(r = 0 ) = 0"*", over a scale that is solely set by 
the confinement radius. We shall note however that the 
numerical simulations underestimate the critical packing 
fraction 0 * at which collective motion occurs, which is 
not really surprising given the number of approximations 
required to establish the interaction parameters in the 
equations of motion Eq. 3. We unambiguously conclude 
from this set of results that Eqs. 1, 2 and 3 include all 
the physical ingredients that chiefly dictate the collective 
dynamics of the colloidal rollers. We now exploit the 
opportunity offered by the numerics to turn on and off 
the four roller-roller interactions one at a time, namely 
the alignment torque. A, the repulsion torque B and 
force 6 , and the dipolar coupling C. Snapshots of the 
resulting particle distributions are reported in Fig. I§. 
None of these four interactions alone yields a coherent 
macroscopic vortex. We stress that when the particles 
solely interact via pairwise-additive alignment torques, 
5 = C = 6 = 0 , the population condenses into a single 
compact polarized swarm. Potential velocity-alignment 
interactions are not sufficient to yield macroscopic 
vortical motion. We evidence in Fig. (top-right and 


bottom-left panels) that the combination of alignment 
(A 7 ^ 0) and of repulsive interactions (5 7 ^ 0 and/or 
6 7 ^ 0 ) is necessary and sufficient to observe sponta¬ 
neously flowing vortices. 

Analytical theory. Having identified the very in¬ 
gredients necessary to account for our observations, we 
can now gain more detailed physical insight into the spa¬ 
tial structure of the vortices by constructing a minimal 
hydrodynamic theory. We start from Eqs. 1, 2 and 3, 
ignoring the C term in Eq. 3. The model can be further 
simplified by inspecting the experimental variations of 
the individual roller velocity with the local packing frac¬ 
tion, see Supplementary Figure 1. The roller speed only 
displays variations of 10 % as 0(r) increases from 10 “^ 
to 4 X 10“^. These minute variations suggest ignoring 
the contributions of the repulsive forces in Eq. 1, and 
solely considering the interplay between the alignment 
and repulsion torques on the orientational dynamics of 
Eq. 2. These simplified equations of motion are coarse¬ 
grained following a conventional kinetic-theory frame¬ 
work reviewed in [4] to establish the equivalent to the 
Navier-Stokes equations for this two-dimensional active 
fluid. In a nutshell, the two observables we need to de¬ 
scribe are the local area fraction 0 and the local mo¬ 
mentum field 011. They are related to the first two an¬ 
gular moments of the one-particle distribution function 
^/(r, p,t) = 7ra‘^{^-6{T — ri)6{p — Pz)), which evolves 
according to a Fokker-Plank equation derived from the 
conservation of ijj and Eqs. and This equation is 
then recast into an infinite hierarchy of equations for the 
angular moments of 'ip. The two first equations of this hi¬ 
erarchy, corresponding to the mass conservation equation 
and to the momentum dynamics, are akin to the continu- 
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ous theory first introduced phenomenologically by Toner 
and Tu [2 E] ; 

dt4’ + foV • ((/in) = 0, (4) 

at((/.n) + uoV-(^(^Q + ^i) =F((/.,n,Q), (5) 

where Q is the usual nematic order parameter. The 
meaning of the first equation is straightforward, while the 
second calls for some clarifications. The divergence term 
on the left-hand side of Eq. is a convective kinematic 
term associated with the self-propulsion of the particles. 
The force field F on the right-hand side would vanish for 
non-interacting particles. Here, at first order in a gradi¬ 
ent expansion, F is given by: 

( 6 ) 

This force field has a clear physical interpretation. The 
first term reflects the damping of the polarization by the 
rotational diffusion of the rollers. The second term, de¬ 
fined by the time rate a = (J^^ 2 a echoes 

the alignment rule at the microscopic level and promotes 
a nonzero local polarization. The third term, involving 
P = (f^^^^r^B(r)dr)/(2a^), is an anisotropic pressure 
reflecting the repulsive interactions between rollers at the 
microscopic level. Eqs.|^and[^are usually complemented 
by a dynamical equation for Q and a closure relation. 
This additional approximation, however, is not needed 
to demonstrate the existence of vortex patterns and to 
rationalize their spatial structure. 

Looking for axisymmetric steady states, it readily fol¬ 
lows from mass conservation, Eq. that the local fields 
must take the simple forms: (/) = 0 (r), 11 = n^(r)e^ 
and Q = Q{r){e^e^ — e^e^), where Q{r) > 0. We also 
infer the relation [—D -1- a(j){l — 2Q)] = 0 from the 

projection of the momentum equation, Eq. on the az¬ 
imuthal direction. This relation tells us that the com¬ 
petition between rotational diffusion and local alignment 
results in a mean-field transition from an isotropic state 
with = 0 to a polarized vortex state with H^ 7 ^ 0 
and Q = | (1 — D/{a(j))). This transition occurs when 0 
exceeds = D/a^ the ratio of the rotational diffusivity 
to the alignment strength at the hydrodynamic level. In 
addition, the projection of Eq. on the radial direction 
sets the spatial structure of the ordered phase: 

2^g-/?(l + 2Q)^=0, (7) 

r dr 

with again Q — | (1 — !<^) in the ordered polar phase. 

This equation has a clear physical meaning and expresses 
the balance between the centrifugal force arising from the 
advection of momentum along a circular trajectory and 
the anisotropic pressure induced by the repulsive inter¬ 
actions between rollers. It has an implicit solution given 
by 


r 


= exp 


2 A-—h A log 



(8) 


0 (r) is therefore parametrized by the dimensionless num¬ 
ber A = (f)'^ ^ reflecting the interplay between self¬ 
propulsion and repulsive interactions. Given the experi¬ 
mental values of the microscopic parameters, A is much 
smaller that unity (A ^ 0.08). An asymptotic anal¬ 
ysis reveals that is the typical core radius of the 
vortex. Eor r < r'^ ^ the density increases slowly as 
0 ^ 0* [1 + for all and Rc. As r reaches 

it increases significantly and then growths logarith¬ 
mically as 0 ^ vq/P \og{r/r'^) away from the vortex core. 
However, r'^ is an integration constant which is solely 
defined via the mass conservation relation: Tri^^^o = 
2'Kr(j){r) dr and therefore only depends on 0 o and R^. 
r* does not provide any intrinsic structural scale, and 
the vortex patterns formed in different confinements are 
predicted to be self-similar in agreement with our exper¬ 
iments and simulations despite the simplification made 
in the model, Eig. [^. In addition, Eq. implies that 
the rollers self-organize by reducing their density at the 
center of the vortex down to (j) = cjA" ^ the mean area frac¬ 
tion at the onset of collective motion, again in excellent 
agreement with our measurements in Eig. li- 

In order to characterize the orientational structure of 
the vortices, an additional closure relation is now re¬ 
quired. The simplest possible choice consists in neglect¬ 
ing correlations of the orientational fluctuations, and 
therefore assuming Q = n^/2. This choice implies that 

n^(r) = (9) 

Eqs. and provide a very nice fit of the experimental 
polarization curve as shown in Eig. Ht, and therefore cap¬ 
ture both the pitchfork bifurcation scenario at the onset 
of collective motion and the saturation of the polariza¬ 
tion at high packing fractions. The best fit is obtained 
for values of and p respectively five and two times 
larger than those deduced from the microscopic param¬ 
eters. Given the number of simplifications needed to es¬ 
tablish both the microscopic and hydrodynamic models, 
the agreement is very convincing. We are then left with 
a hydrodynamic theory with no free fitting parameter, 
which we use to compute the area fraction of the outer 
polarized ring where H^ >0.5. The comparison with the 
experimental data in Eig. [Sfis excellent. 

Eurthermore, Eqs. |^and[^ predict that the rollers are 
on the verge of a phase separation. If the roller fraction 
in the vortex core were smaller (0 < orientational 
order could not be supported and an isotropic bubble 
would nucleate in a polar liquid. This phase separation 
is avoided by the self-regulation of (j){r = 0 ) at 0 *. 


DISCUSSION 

Altogether our theoretical results confirm that the 
vortex patterns stem from the interplay between self¬ 
propulsion, alignment, repulsion and confinement. Self¬ 
propulsion and alignment interactions promote a global 





azimuthal flow. The repulsive iuteractious preveut cou- 
deusatiou of the populatiou ou the geometrical bouud- 
ary aud allow for exteuded vortex patterus. If the rollers 
were uot coufiued, the populatiou would evaporate as 
self-propulsiou iuduces a ceutrifugal force despite the ab- 
seuce of iuertia. 

We close this discussiou by stressiug ou the geuerality 
of this sceuario. Firstly, the vortex patterus do uot rely 
ou the perfect rotatioual symmetry of the bouudaries. As 
illustrated iu Supplementary Figure 2 the same spatial 
orgauizatiou is observed for a variety of couvex polyg- 
oual geometries. However, strougly auisotropic, aud/or 
strougly uou-couvex coufiuemeuts cau yield other self- 
orgauized states such as vortex arrays which we will char¬ 
acterize elsewhere. Secoudly, ueither the uature of the 
repulsive coupliugs uor the symmetry of the iuteractious 
yieldiug collective motiou are crucial, thereby makiug the 
above results relevaut to a much broader class of exper- 
imeutal systems. For iustauce, self-propelled particles 
eudowed with uematic aligumeut rules are expected to 
display the same large-scale pheuomeuology. The exis- 
teuce of a ceutrifugal force does uot rely ou the directiou 
of the iudividual trajectories. Shakeu rods, couceutrated 
suspeusious of bacteria, or motile biofilameuts, amoug 
other possible realizatious, are expected to have a simi¬ 
lar phase behavior. Quautitative local aualysis of their 
spatial patterus [3 EHH US] would make it possible to 
further test aud elaborate our uuderstaudiug of the struc¬ 
ture of coufiued active matter. For iustauce, the polar 
order fouud iu coufiued bacteria is destroyed upou iu- 
creasiug the size of the coufiuemeut. The aualysis of the 
spacial distributiou of the bacteria could be used to gaiu 
iusight ou the symmetries aud the maguitude of the ad- 
ditioual iuteractious mediated by the host fluid which are 
respousible for the emergeuce of bacterial turbuleuce [16] . 

Iu couclusiou, we take advautage of a model exper- 
imeutal system where eusembles of self-propelled col¬ 
loids with well-established iuteractious self-orgauize iuto 
macrosopic vortices wheu coufiued by circular geomet¬ 
ric bouudaries. We ideutify the physical mechauism that 
chiefly dictates this emergeut behavior. Thauks to a com- 
biuatiou of uumerical simulatious aud aualytical theory, 
we demoustrate that orieutatioual coupliugs aloue cauuot 
accouut for collective circular motiou. Repulsiou betweeu 
the motile iudividuals is uecessary to balauce the ceu¬ 
trifugal flow iutriusic to auy ordered active fluid aud to 
stabilize heterogeueous yet mouophasic states iu a broad 
class of active fluids. A uatural challeuge is to exteud 
this descriptiou to the compact vortices observed iu the 
wild, e. g. iu shoals of fish. Iu the abseuce of coufiu- 


iug bouudaries, the ceutrifugal force has to be balauced 
by additioual deusity-regulatiou mechauisms [35] |36] . A 
structural iuvestigatiou akiu to the oue iutroduced here 
for roller vortices could be a powerful tool to shed light 
ou deusity regulatiou iu uatural flocks which remaius to 
be elucidated. 

METHODS 

Experiments. We use fluoresceut PMMA colloids 
(Thermo scieutific G0500, 2.4/rm radius), dispersed iu 
a 0.15 mol.L“^ AOT/hexadecaue solutiou. The suspeu- 
siou is iujected iu a wide microfluidic chamber made of 
double-sided scotch tapes. The tape is saudwiched be¬ 
tweeu two ITO-coated glass slides (Solems, ITOSOL30, 
80 um thick). Au additioual layer of scotch tape iu- 
cludiug a hole haviug the desired coufiuemeut geometry 
is added to the upper ITO-coated slide. The holes are 
made with a precisiou plottiug cutter (Graphtec robo GE 
6000). The gap betweeu the two ITO electrodes is cou- 
staut over the eutire chamber H = 220 /im. The electric 
field is applied by meaus of a voltage amplifier (Stauford 
Research Systems, PS350/5000V-25W). All the measure- 
meuts were performed 5 miuutes after the begiuuiug of 
the rolliug motiou, wheu a steady state was reached for 
all the observables. 

The colloids are observed with a 4X microscope ob¬ 
jective for particle trackiug, PIV aud uumber-deusity 
measuremeuts. High-speed movies are recorded with a 
GMOS camera (Easier ACE) at a frame rate of 190fps. 
All images are 2000x2000 8-bit pictures. The particles 
are detected to sub-pixel accuracy aud the particle tra¬ 
jectories are recoustructed usiug a MATLAB versiou of 
a couveutioual trackiug code The PIV aualysis was 
performed with the mpiv MATLAB code. A block size 
of 44 jam was used. 

Numerical simulations. The simulatious are per¬ 
formed by uumerically iutegratiug the equatious of mo¬ 
tiou Eqs. 1 aud 2 . Particle positious aud rolliug direc- 
tious are iuitialized raudomly iuside a circular domaiu. 
lutegratiou is doue usiug au Euler scheme with au adap¬ 
tive time step St, aud the diffusive term iu the equatiou 
for the rotatioual dyuamics is modeled as a Gaussiau 
variable with zero meau aud with variauce 2D/St. Steric 
exclusiou betweeu particles is captured by correctiug par¬ 
ticle positious after each time step so as to preveut over¬ 
laps. Bouuciug off of particles at the coufiuiug bouudary 
is captured usiug a pheuomeuological torque that reori- 
euts the particles towards the ceuter of the disc; the form 
of the torque was choseu so at the reproduce the bouuc¬ 
iug trajectories observed iu the experimeuts. 
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Supplementary Information 

Supplementary Figures 



FIG. 5. Variation of the roller velocity with the local area fraction. EjEQ = 1.4. As (j) varies from 10 ^ to 4 x 10 vq only 
increases by ~ 10%. 






1mm 

FIG. 6. Superimposed pictures of colloidal rollers forming a vortex patterns in polygonal geometries. Packing fraction: 
00 = 3 X 10~^. EjEQ = 1.4. The self organization of the population is robust to the shape of the confinement. These pictures 
show that the same heterogeneous vortex pattern emerge in polygonal confinements which do not have a perfect rotational 
symmetry. The center of the vortex is dilute and less ordered that the outer region where polar order is pronounced. The shape 
of the polygonal geometries is reflected by the velocity and density prohles yet the salient features of the vortex is preserved. 
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xlO^ 



FIG. 7. Mean-squared displacement of the rollers, {[ri{t + T) — plotted as a function of the lag time T. Open symbol: 

experiments (error bars ISD). Full line: theoretical prediction. 


Supplementary Discussions 

Supplementary Discussion 1 | Dynamics of a single roller 

At low packing fraction, 00^0*? Ihe colloids behave as non-interacting persistent random walkers. Their motion 
is described by: 


dtVi = i;oPi = vo{cos0i,sm0i), (10) 

dt0i=Ci{t). ( 11 ) 


where is a white noise with zero mean and variance = 2D5(t — t'). The velocity autocorrelation function 

decays exponentially as 


{voVi{t + T) ■ vo'Pi{t)) = Vq e 


( 12 ) 


and the mean-squared displacement is given by: 

+ T) - r,(t)]2) =2^{DT-1 + e-^^) . (13) 

The later expressions are used to fit the experimental data, Fig. 2E (main text) and Supplementary Figure 3 below, 
and provide the following values for the particle speed and diffusivity; 


no = 493 ± 17 /rm.s 
D-^ = 0.31 ± 0.02 s. 


(14) 

(15) 
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Supplementary Discussion 2 | Population of interacting rollers 

Microscopic model. In [T], we have theoretically described the microscopic dynamics of a population of colloids 
rolling on a conducting surface. We briefly summarize this model. A colloid powered by the Quincke mechanism, 
when rotating close to a wall, exchanges momentum with this solid surface and translates at a speed given by 


= 


lirT 



- 1 . 


(16) 


Eq is the critical electric field below which the particle does not rotate. The Maxwell-Wagner time r characterizes 
the dynamics of the electric charges at the colloid surface, which are responsible for the Quincke instability, and 
jit are mobility coefficients accounting for the viscous drag exerted by the liquid in the vicinity of solid wall that 
depend logarithmically on the distance to the solid wall [T]. In order to account for the roller-roller interactions 
at long distances, we have computed the electrostatic and hydrodynamic fields induced by the motion of a colloid. 
Assuming pairwise additive interactions, we have shown that these fields promote the alignement of the velocity of 
neighboring particles. Within this framework, the speed of the colloids is taken to be a constant as it relaxes to vq 
much faster than the typical timescale of the orientation dynamics. In addition to these far-field couplings, we model 
the short-distance repulsion between colloids by an effective hard-core potential. The resulting equations of motion 
are given by Eqs. [I]-[3] in the main text, where the functions A(r), 5(r), C{r) are given by: 


A(r) = Ai 1 

(°) e(’') + ^2( 


B{r) = Bi 1 



C{r) = Cl 

to 

to 

+ 

(^)’e(.) 


+ c.0'e(r). 


(17) 

(18) 
(19) 


These electrostatic and hydrodynamic couplings are exponentially screened over a distance set by the channel height 
H (see Fig. I A, main text). For sake of simplicity, we approximate the screening function 0(r) by the step function 
0(r) = I if r < H/tt and 0(r) = 0 otherwise. In addition, the coefficients of the above functional forms are given by: 


Ai 

A 2 

Bi 

Cl 

C 2 






6r-q^-l 

/^r 


Ai, 


1^ 


El)' 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


fir and /ix are mobility coefficients which only depend on the viscosity of the liquid and the gap d between the 
colloid and the surface, depends on the dielectric permittivities ei and Cp of the liquid and of the particles: 
= {ep — € 1 ) {€p + 2ei) The derivation of this model is provided in the supplementary informations of [1]. 


Estimation of the simulation parameters. In order to perform numerical simulations relevant to our 
experimental conditions, we have estimated the coefficients of the equations of motion as follows. The speed vq 
and the rotational diffusivity D have been deduced from the single-particle dynamics, Eqs. ([T^ and ([T^. The 
threshold electric field Eq is measured experimentally as the critical value at which the colloids start moving. We 
use typical values for the dielectric permittivities of hexadecane (e/ = 2.2eo) and PMMA colloids (e^ = 2.6eo) [2] 
to evaluate = 0.06. The mobility coefficients are estimated. We assume the distance between a particle and 
the surface to be d ^ 50 nm. Although this parameter is not controlled precisely, it only yields small corrections 
to the mobility coefficients in the limit d ^ a, and weakly impacts the particle dynamics. Using the expressions 
derived in nn, we find jig = 0.30, ii^/iir — 1-b and jltl[ir — 8.7 x 10 Finally, the Maxwell-Wagner time r was 
calculated from Eq. (16): r = 0.29 ms. As a result, we obtain the following values for the microscopic coefficients: 
Ai = Bi = Cl = 0.9 r^, A 2 = I.0r“^, and C 2 = 1.7 









14 


[1] Bricard A., Caussin J.-B., Desreumaux N., Dauchot O. & Bartolo D. Emergence of macroscopic directed motion in popu¬ 
lations of motile colloids. Nature 503, 95-98 (2013). 

[2] Pannacci N., Lemaire E. & Lobry L. DC conductivity of a suspension of insulating particles with internal rotation. Eur. 
Phys. J. E28, 411-417 (2009). 

[3] Goldman A. J., Cox R. G. & Brenner H. Slow viscous motion of a sphere parallel to a plane wall - I Motion through a 
quiescent fluid. Chem. Eng. Sci. 22, 637-651 (1967). 

[4] Goldman A. J. , Cox R. G. & Brenner H. Slow viscous motion of a sphere parallel to a plane wall - II Couette flow. Chem. 
Eng. Sei. 22, 653-660 (1967). 

[5] O’Neill M. E. & Stewartson K. On the slow motion of a sphere parallel to a nearby plane wall J. Eluid Meeh. 27, 705-724 
(2006). 

[6] Liu Q. & Prosperetti A. Wall effects on a rotating sphere J. Eluid Meeh. 657, 1-21 (2010). 

[7] Desreumaux N., Caussin J.-B., Jeanneret R., Lauga E. & Bartolo D. Hydrodynamic fluctuations in conhned emulsions. 
Phys. Rev. Lett. Ill, 118301 (2013). 



